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Abstract 

] We find the density matrix corresponding to the vacuum state of a massless Dirac field in two 

■ dimensions reduced to a region of the space formed by several disjoint intervals. We calculate explicitly 
its spectral decomposition. The imaginary power of the density matrix is a unitary operator imple- 
menting an internal time flow (the modular flow). We show that in the case of more than one interval 

O || this evolution is non-local, producing both, advance in the causal structure and "teleportation" be- 

tween the disjoint intervals. However, it only mixes the fields on a finite number of trajectories, one for 
each interval. As an application of these results we compute the entanglement entropy for the massive 
multi-interval case in the small mass limit. 

1 Introduction 

An observer with constant self-acceleration a in vacuum Minkowski space-time sees himself immersed 
in a thermal bath at temperature T = a/(2ir) [lj. This is the Unruh effect. One explanation for this 
phenomenon is that the constant acceleration impedes this observer to be in contact with all the Minkowski 
space. In fact, there is a portion of the space-time (a Rindler Wedge), including half of a spatial hyperplane, 
from which the signals cannot reach the observer. In consequence, the degrees of freedom lying in this 
region are not relevant for the description of the physics in the observer's world. If we trace over them, 
the vacuum state becomes a mixed state on the observer's half space V. This mixing reflects the presence 

■ of quantum entanglement between the seen and the hidden regions. 

. In a more general case, we can consider localization in a different region V. The corresponding local 

state py is defined by the property that it leads to the same expectation values for the operators localized 
in V than the vacuum state. In presence of a regularization, py can be written in terms of a density 
matrix, in such a way we have 

tr(p v O) = (0\O\0) , (1) 



ON 

o 



X 



for any operator O localized in the region V (or more precisely localized in the domain of dependence of 
V). In this paper we will assume implicitly such a regularization (a lattice regularization for example), 
which will allow us to use the language of local density matrices (in the continuum limit the local algebras 
of operators are a special kind of von Neumann algebras, and the states on the local algebras should be 
treated with more sofisticated mathematical machinery than density matrices). 
Any density matrix can be written as an exponential 

pv = ce- n , (2) 

of a Hermitian operator TC, and where c is a normalization constant. The thermal interpretation of the 
local state py arises when one considers an internal time r, with unitary evolution given by the operator 
Py, inside the algebra of the fields in the region V. With respect to this notion of time evolution, py is 
indeed thermal, with temperature rescaled to 1. 
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In the algebraic axiomatic approach to quantum field theory, the internal time evolution is called the 
modular flow and the Hermitian operator corresponding to 7i the modular Hamiltonian [2J. These are 
well defined using the Tomita-Takesaki theory □, which plays an important role in the structure. The 
ciclic and separating vector state of the Tomita-Takesaki theory (here the vacuum state |0)) is related to 
pV through eq. ([I]). 

In a different context, this mechanism of time flow generation from a thermodynamical state has 
been related with the time problem in quantum gravity [I]. The thermal time hypotesis aims to explain 
the emergence of a physical time flow in a general covariant quantum theory as determined by the 
thermodynamical state rather than been a universal property of the mechanical theory. 

In the Unruh effect the set V corresponds to a Rindler wedge. This particular case has a remarkable 
feature. The Hamiltonian has the universal expression H = 2ttK, in terms of the boost operator K which 
keeps the Rindler wedge fixed. This is valid for any quantum field theory [5]. The internal evolution is local 
and causal. It is given by point transformations along the integral curves of the boost operator, which are 
precisely the trajectories of constant acceleration observers. In curved space-time a similar universality 
occurs for the black hole, where the local hamiltonian for the outside region and the Hartle-Hawking state 
is proportional to the time translation operator for asymptotic observers. 

For different type of sets these features should be lost in general. An exception is the internal evolution 
corresponding to spherical sets V (which are single intervals in two dimensions) in conformal theories. 
This is also local, and given by the one parameter group of global conformal transformations which keep 
V fixed 0E]. However, no other explicit example of internal evolution is known in quantum field theory. 
In particular no non local example has been explicitly computed. 

In this work we provide one such example. We calculate py when V is formed by several disjoint spatial 
intervals in two dimensional Minkowski space, for a massless Dirac field. We explicitly diagonalize the 
reduced density matrix in this case, and obtain the internal evolution. We find some surprising features. 
Even if the internal evolution is non local, it does not mix all the operators on the Cauchy surface, but the 
operators in selected points on each interval. Also, we find that the total posible teleportation between 
the different regions appears to be limited. 

In curved space a similar phenomena should be present in the case of several black hole regions. In 
this context, the present example provides an exactly solvable toy model in order to explore the subtle 
effects of the multipartite entanglement present on the vacuum fluctuations, and their relation with the 
geometry. Also, these effects should give place in flat space to correlations on the meassurements of two 
or more accelerated observers. 

As an application of the results, we also calculate the multicomponent entanglement entropy in the 
massive case under some approximations. We show the property of extensivity of the mutual information 
which holds for the massless case is lost in the massive one. 

2 Structure of the density matrix 

An expression for the reduced density matrix py in terms of the two point correlators is known for free 
Dirac fields E]. This follows from the fact that the expectation values of polynomials on the fields 
located in V computed with the help of py must obey Wick's theorem (here the global state is the 
vacuum). This fixes py to be of the form ^ where TL is a Hamiltonian for V which is quadratic on 
the fields inside this region. For discrete bosonic systems analogous expressions for TL exist which are 
quadratic on the field and momentum operators. However, the generalization of these expressions to the 
continuum case is not known |10j . 

In this paper we are interested in free Dirac fields in two dimensions. However, the formulas in this 
section apply to any dimensions with the obvious modifications. Let us consider the general case of regions 
which are arbitrary spatial curves, not necessarily included in a single spatial line, as shown in figure 1. 

1 In the formalism of quantum fields, the Tomita-Takesaki theory can be interpreted in terms of thermofield dynamics [3]. 
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Figure 1: Two spatial sets V and V' in two dimensional Minkowski space. They have the same causal 
domain of dependence which is the diamond shaped set in the figure. Light rays are shown at ±45°. 



Thus we write H as 

H= / dsi ds2^\{si)H ij {si,s 2 )^j{s2) ; #?j(s 2 , s x ) = -Hy(si, s 2 ) , (3) 
Jv 

where ^(s) = ^>(x^(s)), and s%, s 2 are length parameters along the spatial curve V. The Dirac field ^f(x) 
satisfies the Dirac equation 

(iffy-m)* = 0, (4) 
and the canonical (equal-time) anticommutation relations 

{^ i (x,t),^ j (y,t')]\ t=ti =8(x-y)6 ij . (5) 

For the fields ^f(s) parametrized with the length s on the surface V the slope of the curve (its relative 
Lorentz boost with respect to the t =constant lines) has to be taken into account, and equation © 
generalizes to 

{*i(si), ^j(s 2 )} = 5(ai - s 2 )5ij , (6) 

where we have written 

%( S ) = ^(s)fj(s) , fj(s) = 7 °7%(s) , (7) 

and rj^(s) the unit vector normal to the surface V at the point s. Note that fj(s) = fj(sy is hermitian and 
fj(s) 2 = l. 

Write the two point correlation function restricted to V as 

C ij -(ai,«2) = <0|*i(ai)^(a 2 ) t |0> • (8) 

Then we must have 

tr (py^(si) tft( a2 )) = djisi, s 2 ). (9) 

Taking into account the expressions ([2]) and (|3j) for the density matrix and the field algebra © on 
the surface, the equation Q leads to the following relation between the Hamiltonian kernel H and the 
correlator 

H(si,s 2 ) = -fj(si) log [{Cf))- 1 - 1)] (si,s 2 ) . (10) 
One can also summarize this in a more compact manner, 

p v = ce~ H , K= ds 1 ds 2 ^(s 1 )H(s 1 ,s 2 )^(s 2 ) , (11) 

Jv 

with 

H = - logiC- 1 - 1) (12) 

and 

C(ai,s 2 ) = (0|*(fl 1 )#(a 2 )|0) . (13) 
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2.1 Causality 

The explicit form of the density matrix in terms of the operators lying on any spatial curve allows us to 
show that the reduced density matrix does not change with the Cauchy surfaces having the same domain 
of dependence (see figure 1). Heuristically, this means the information cannot travel faster than the light 
speed, and in consequence it must be the same in the surfaces V and V in figure 1. In order to show 
this, consider the equation 

*i(z) = J ds fax), ^(s)} ^(s) , (14) 

which is valid for any space-time point x in the domain of dependence of V, and where s is again a distance 
parameter along the curve V. This follows by the uniqueness of the solution of the Dirac equation since 
the right hand side satisfies the Dirac equation in x, and on the surface V the equation is satisfied due to 
([6]). Then, let us call 

G^( Sl , S2 ) = {^( Sl ),^( S2 )} , (15) 

where ^i(si) is a field on V' and ^fj(s2) is a field on V. We have the following identities in operator 
notation 



where 

G vv ' =fj v (G v ' v yfj v ' . (19) 
Now, considering the anticommutator in the surface V 

iT = Uv^y'X = G v ' v h? v ,$ v } G vv ' = G V ' V G VV ' . (20) 

Thus, G v ' v and G vv ' are inverse to each other. It is then easy to see that the similarity transformation 
of the correlator (|16M18p cannot change Ti in eqs. (|imi3p . This means that we can use any spatial slice 
for a given causally complete set (that is, a diamond shaped set arising as a domain of dependence of a 
spatial set) in order to compute the density matrix. 



3 The massless case 



®V = V V G VV ', 
WV _ r<V'V nVV riW 



(16) 
(17) 
(18) 



In this paper we focus on a Dirac field in 1 + 1 dimensions. The correlator for two points x and y with 
spatial separation [x — y) 2 < is 

C(x,y) = ]d, [9(x-y) 2 S gn(x°-y )} fV + ^K (m\x - y\) 7 ° - ^ Kl (m\x - y\) { * ~ % 7 ^ 7 ° , (21) 
4 L J 2n 2tt \x — y\ 

where K a (x) is the standard modified Bessel function. The first term gives just l/26(s\ — s 2 ) in C for 
any surface. 

In the massless case the problem can be decomposed into the different chiralities. Let us define the 
chirality operator 7 3 = 7°7 1 , with 7 3 = ( 7 3 )^ and ( 7 3 ) 2 = 1, and the chirality projectors Q± = (l± 7 3 )/2. 
Let us also introduce the null coordinates 

u±=t±x. (22) 

The massless Dirac equation leads to d u± ^^z = 0, where *$>± = Qzp^f are the field components of definite 
chirality (note that ^± has chirality ^1). Then we have 

* = *+(n + ) + ^_(n_). (23) 
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These fields depend only on one of the null coordinates, and satisfy the anticommutation relations 

{*±(4),vl4(4)} = diui-ul), (24) 
{*±(4),4(4)} = 0. (25) 

Using this decomposition and the massless limit of the correlator (I21D in eqs. (|iml"3l) we have (for any 
choice of surface) 

p = p(u+) ®p(w_) = c e~ n + (g>e' H - , (26) 
H = H++H-, (27) 

H± = f duidi4*i(u±)-ff±(t4>^)*±(«±)> ( 28 ) 
fl± = - logpi 1 - 1) . (29) 

Here V± are the projections of V in the u + and u_ coordinate axis respectively, and D± is now the scalar 
kernel 

DfevJ.^-rt-JL^i-j (30) 

where the domain of the real variables cc,y is V±. Here and in the following the distribution (x — y)~ is 
taken with the principal value regularization. 

3.1 Eigenvectors and the resolvent of D 

Therefore, in the massless case, the problem of finding an explicit expression for the local density matrix 
and related quantities is equivalent to the one of solving the operator with kernel D. Fortunately, the 
resolvent of the kernel (x — y)~ l in multicomponent subsets of the real line is known from the theory 
of singular integral equations of the Cauchy type [11]. We write in a generic way a set V = (ai,b\) U 
(02, 62) U ... U (a n , b n ), formed by n disjoint intervals, where a, < 6j, b{ < dj+i. For a set of intervals Ii in 
spacetime, ordered from left to right in the spatial coordinate, we write the («+,«_) coordinates of the 
left and right extreme points of the intervals as (ujf^u^) and (u+^u^). Thus in the following (aj,6j) 
may mean either (af,bf) = (n^,M^) or (a^,b~) = (u R ^ n+1 _^, u^ n+1 _^), according to where we are 
referring to the u + or the n_ coordinate axis (note that in this later the order of the points is inverted) . 
In this section we use generically the variables x, y as meaning either coordinates on the null axis or 
in the axis u_. 

The resolvent of D, 

R°(P) = (D — 1/2 + /3)" 1 = (-^ — +p S (x- y)) 1 , (31) 

V 2vr x - y ) 

has the following expression [11 J 

, / ,• p -^ lo s(fe^)(-(-)-^))\ 
R\(3) = {(3 2 - 1/4) hs(x - y) + — — , (32) 

where 

^M-m^i)' (33) 

Moreover, we have been able to obtain the complete spectrum of (x — y)~ l in V. We have guessed 
the solutions for the eigenvectors (somehow inspired in some formulae of the book [llj, and also using 
numerical integration) and checked directly they where the right eigenvectors, by integration. 
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It has n degenerate eigenvectors 

I dy— **(y) = A.^( l/ ), (34) 
k = 1, ..,re, corresponding to the eigenvalue 

A s = i7rtanh(7rs) (35) 



for any s S R. The eigenvalue corresponding to H is 2tts. Eq. (|35p and most of the identities below can 
be proved by writing the integrals in the complex z plane. 

The eigenvectors can be chosen to form an orthonormal basis as 



V k (x) = XX ^ fc -—^ e ~ tsz ^ , (36) 

N k (Xl^lx-aiWx-bi" 1/2 



where 6v(x) is a function which takes the value (— 1) J+1 in the j th interval (aj,bj). The normalization 
factor is 



" " ^ I ^ (6i " «*) IL^fo " 6i) IE=i(a* " &i) 



2?r 2tt 

i=l 



fc=i 

We also have the useful relations 



and 



3.2 Modular Hamiltonian 

The Hamiltonian kernel H = — log(D _1 — 1) writes in terms of the resolvent (|31j) as 



Changing the variable (3 to a = ^ log( o.M^ )) we have from (f32l) 



1 -*(*(*)-*(»)) = _ 2m 6(z(x)-z(y)) 



(37) 



We have from this 

fx 2 ^"-Sr (38) 



where Lt is the sum of the n interval lengths. The vectors (J36]) then satisfy 

/ dx* k s *(x)y k s :(x)=5 kjk ,6(s-s'), (39) 
Jv 

and form a complete orthonormal basis 

n POO 

ds^ k s *(x)^ k s (y) = S(x-y). (40) 

7. -I J-OO 



/ cfe^(z) = (-l)" +1 7rsech(vrs)iV fe - 1 , (41) 
Jv 

Im / dxx^ k {x) = (-l) n svrsech(vrs)iV fe - 1 Li . (42) 



£T = - / dp (H?(J3) + B?(-0)) . (43) 

Jl/2 



H = -i da e -vr(*W-z(v» = _ 27ri — Z^Zl (44) 

J-oo x-y x-y 
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Figure 2: Three related trajectories for a three interval case. The null coordinates in the three diamonds 
are u±(z±), I = 1, 2, 3, and z± 6 (— oo,oo), with z + — Z- =constant. The related points do not show any- 
geometric symmetry in general. 

The expression in (|44j) has two different contributions, 

H = H\ oc + H no \ oc . (45) 

H no i oc is non local, and occurs when z(x) = z{y) but x ^ y. The local one H\ oc is due to the solution 
x = y of the equation z(x) = z(y). After changing the variables in the argument of the 5 function, this 
last contribution gives a product of distributions of the form 8{x — y)/{x—y). The different regularizations 
of this quantity must vanish outside x = y, and give 5(x — y) when multiplied by the function (x — y). 
The only solutions are of the form —6'(x — y) + kS(x — y), where k is regularization dependent. Here, k 
is fixed by requiring that H is hermitian. This gives for the local term 

The non local term follows by solving the equation z(x) = z{y) for x ^ y. The function z(x) increases 
monotonically inside each of the intervals Ij, being — oo on the points a, and tending to oo on the other 
end point 6j of the interval. Thus, it follows that given any z £ (—00,00), there is exactly one solution 
xi(z) € of (]33p for each interval. These solutions are given by the roots of a polynomial. Note that the 
different points xi(z) which are related to each other tend to the left (right) of the respective intervals 
simultaneously. In the case of two intervals the points connected to each other in the Hamiltonian are 

, x _ (hh ~ Qiag) x 2 {z) + (61 + 6 2 )aia 2 - (ai + 02)^2 (A7 , 

\ Z J til \ f \ IT * \ ) 

(01 + o 2 - ai - a 2 ) X2{z) + aia 2 - 0i& 2 

This is a global conformal transformation of square 1, sending a\ to a 2 and b\ to 6 2 . 
The non local part of the Hamiltonian then writes 

1 ( d \ — ^ 

H noloc = -2m Yl 7^)(x S(y- Xl (z(x))). (48) 

Remarkably, this non local Hamiltonian only mixes a finite number of points, the Xi(z), one from each 
interval, having the same z. 

It is possible to obtain the same expressions for H using directly the decomposition in terms of the 
eigenvectors 

/oo 
dss* k s (y)y k s *(x). (49) 
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3.3 Internal dynamics 

The internal dynamics inside V (one parameter group of modular transformations) is given by the unitary 
operators p lT . In order to see how these transformations act we look at their action on the field operators 

V(x,t) = p iT t>(x)p- iT . (50) 

Using the expressions (|26H29|) of the density matrix in terms of a quadratic Hamiltonian and the anticom- 
mutation relations (|24H25[) . we have that the chiral components ^±(u±,t) satisfy a Dirac-type equation 
with Hamiltonian H±, 

-jf = -i[H±,* ± ]=iH±*±. (51) 

For one interval it is known that H = H\ oc is proportional to the generator of global conformal 
transformations which keep the interval fixed. Eq. (|5ip leads in this case to the point transformation for 
the operators [21I51IT2] 



-M« ± (t),t) = J &iz^)fezj4) , ±( 4), (52) 

V ( U ±( T ) - a±)(b± - u±(t)) 
, s _ b ± (ul-a ± ) + e- 2 ^a±(b±-ul) 
U±[T) ~ (ul-a±)+e-^(b±-ul) ' { } 

where (a±, b±) are the projections of the interval in the u± coordinate axis (a+ and 6_ correspond to the 
left endpoint while a_ and 6+ correspond to the right one). 

For more than one interval the solution of (|51p can be given by expanding in eigenvectors of H±. 
However, these expressions are rather obscure. In order to understand the local dynamics let us write 
the field in the related points in the different intervals as a vector of fields, ^i(z) = ^>{xi{z)). Then 
eq. (|5ip writes (the same equation applies for both chiralities, we omit the ± subindices for notational 
convenience) 

dr dz dz \ dz J 1 dz xi — xy 

The characteristic trajectories z(t) of the first order differential equation are 

z(t) = Z0 + 27TT, (55) 

for a constant zq. Thus, the fields ^/(r) = ^ / ;(z(r),r), will mix to each other for different / along the 
characteristic. Only the n trajectories xi(z(t)), with a given fixed zq will mix to each other, and the 
movement is such that z moves at constant speed 2-7T with respect to r. The figure (2) shows a group of 
three related trajectories inside the domain of dependence of three intervals (the three diamond shaped 
sets). In figure (3) the flux of trajectories are shown for two diamonds in three different relative positions. 

On a group of related trajectories the ^(r) evolve with r according to a coupled set of ordinary 
differential equations. Writing 

^(r) = (<fc,/(fe) 1/2 tfi(T), (56) 

we have 

'dxj\2 ( dx v \ 2 



^ = 2,Vi^ \±L 9 V , (57) 
ar X] — X]> 

where X\ and z are functions of r according to (|55p . Note that since the matrix on the right hand side is 
antisymmetric, the sum ^2i(^i(t)) 2 is a constant, and ^(t) = O(r)^ r (0), with 0{t) orthonormal. 
The analytic solution of this equation for n = 2 can be obtained explicitly. In this case we write 

V sm(fc') cos(y) I 
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Figure 3: Two intervals of size d separated by: (a) (l/4)cZ (top pair of diamonds), (b) (l/70)d, and (c) 
(1/1000)g? (bottom pair of diamonds). The trajectories in (a) approach the ones corresponding to single 
independent diamonds. In (b) the curves get distorted by the proximity of the diamonds. In (c) the 
curves are very similar to the ones of a single diamond of size 2d (shown with dashed lines), excepting 
when they hit the null boundaries and run along them afterward. The trajectories shown for the two 
diamonds are equal spaced in the t = plane, and they are not related to each other. 



and the solution of the differential equation (|57p gives 

tan(fl) = (-i-4)(bi + b 2 -a 1 -a 2 ) (gg) 
V ( & i - oi)(«2 - &l)(&2 - oi)(&2 - a 2 ) 

The tangent of the rotation angle increases linearly with the coordinate x\ as it moves from x\ to b\. 
It does not diverge, meaning that the rotation of components never completes, there is no complete 
"teleportation" or swap between the components. 



4 Entropy 

From (|35h it follows that the spectrum of py is continuous and covers all the interval (0, 1). This means 
that even if py is well defined as a state on the operator algebra, it is not a bona-fide density matrix, 
since it is not trace class. Then, the expression for the entropy which follows from (|llH13p . 



S(V) 



-tr 



;i-(7)log(l-(7) + Clog(<7) 



(60) 



has to be regularized. We will find more convenient to express the entropy in terms of the resolvent 
R = (C - 1/2 + (3)- 1 as 



S(V) 



d(3tr 



1/2 



(0-1/2) (R(P)-R(-P)) 



2/3 



0+1/2 



(61) 
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With the knowledge of the spectrum and the resolvent for the massless case we are in position to compute 
the entropy and to do some expansions for the massive case. 



4.1 Massless case 

In this case we use directly the chiral decomposition and the resolvent (I32p . In applying (I6ip the term 
proportional to the identity cancels with the corresponding one in the resolvent. Then we have, for any 
chiral component, a contribution of the form 



1 



oo i- sin 



S(V) = -- I d(3 I dx lim l2w ; , (62) 

1 ' Wi/2 Jv W + 1/2) (x - y) v ; 

where we are using the notations of sections 3.1 and 3.2, and here x, y should be understood as coordinates 
on one of the chiral axis. One can integrate first over f3, giving 

S (v) = r * run c ° th <<f 1 - W 1 = 1 = - / *> t (— - — 

Jv y^ x {x - y){z{x) - z{y)) 12 J v f-£\x-ai x- 

= g I ^log|6j - ai\ - ^log|aj - aj\ - ^log|6j- bj\ -nloge J . (63) 

\ ij i<j i<j j 

Here e is a distance cutoff introduced in the last integration, which extends from ai + e to hi — e in 
each interval. The divergent terms on the entanglement entropy are associated to the boundary of the 
region, where the high frequency modes give place to an unbounded entanglement with the outside 
region. Physically sensible cutoff independent quantities can be obtained from the regularized entropy by 
appropriate subtractions. For example, the mutual information I(A, B) = S{A) + S{B) — S{A U B) is 
finite and well defined, since the boundary terms cancel in this combination. 

This same result was obtained by bosonization in [13]. The formula (|63p applies for each chirality. 
The same formula but with coefficient 1/3 instead of 1/6 gives the entropy for a set V on a spatial line 
(two chiralities). 

4.2 Massive case: small mass limit 

For the massive case we consider only sets lying in the same spacial line, and thus the correlator reads 
C(x,y) = -5{x - y) + —K (m\x - y\)j° + — #i(m(x - y)) 7 3 . (64) 

The short distance expansion of the correlator up to second order in the field mass is C = Co + C\ + 
Ci + ... where 

C (x,y) = l5(x-2/)JI + ^-—L- 7 3 , (65) 
Cx(x,y) = -^(lE + logi^^-))^, (66) 



2vr V V 2 

C 2 (x,y) = ^- y )f 7E -i+logf^^))- ! . (07) 



4?r v V 2 

and "fE is the Euler constant. The perturbative expansion of the resolvent for small mass is done by 
considering the formula 

R V (f3) = Ry(f3) - RyiftdRy^) - Ry{f3)C 2 Ry(f3) + RyifrCiRyifrdRyiP) ~ ... . (68) 
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The first term gives the massless contribution (I63p . The second one is of first order in the mass, and does 
not contribute to the entropy since it is traceless. From (|65H68p the expansion of the entropy has the form 



oo 2i 



s = J2Yl Si >i m2i logj n • ( 69 ) 

Let us calculate the contributions which are second order in the mass coming from the third and 
fourth terms. First we calculate the leading log 2 (to) term. We find convenient to use the expression of 
the resolvent in terms of the eigenvectors ^f k (x) of (x — y)~ l ■ From ([65]) we have 

n r+oo 

R° v (f3)(x,y) = Y^ / d8* k 8 (x)M(0,s)* k a *(y), (70) 

where 

M(P,s)=^I- t -^pl^y 1 . (71) 
Thus, the log 2 (to) contribution which comes from the fourth term in (1680 is 

2 poo poo poo 

6S 2 ,2 = S2,2^ 2 log 2 (TO) = -■— ^ bg 2 (TO) / d(3 ds d S ' {(3 - 1/2) 

\ / J 1/2 J — 00 J — 00 

N N 

£ £ |A(a, k)\ 2 \A(s>, k')\ 2 tr(M(/3, s) 7 °M(/?, s') 2 7°) , (72) 



k=i k'=i 

where 

A(*,fc) = / ds**(a;). (73) 
Jv 

From (J3SJ) and flU]) we have 

AT 

^|A( S ,A;)| 2 = ^sech 2 (^)Li. (74) 



fc=i 

Using this, the integrals in (1721) can be found exactly giving 

<^2,2 = -!^log 2 (TO). (75) 

The next contributions are of order to 2 log (to). They come from the third and fourth terms in (|68|) . 
We can divide this contribution in three terms, 

SS 2 ,i = s 2 ,i to 2 log(m) = Ai + A 2 + A 3 . (76) 

The contribution coming from the terms in C\ which does not contain log \x — y\ we call Ai, the one from 
C 2 , we call A 2 , and the last one A3 which involves log \x — y\ in C\. Ai is readily evaluated since the 
relevant calculation is the same as above, 

a lo §( 2 ) - HE 2 T 2 , / \ 

Ai = to L t log(m) . (77) 

The contribution of C 2 to this order is 

,2 



2tt 



log(m) / d/3 / ds(/3- l/2)tr( 7 d M(/3,s)^) 



1/2 J-00 

,\ 

x 

fe=i 



Y, [ dx [ dy{x-y)* k a *(x)V k a (y)- (78) 
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Using (jSHD, (35) and <$2§ we find 
Finally, the expression for A3 is 



A 2 = -- L t 2 m 2 log(m) . (79) 



TV N 



z poo poo poo 

A 3 = -—z log(m) / d/3 / ds / ds' V V (J3 - 1/2)A(«, k)A(s', k')* 

(tr(M(/3, s ) 2 7 ^(/5^')7°)-(/3 /dx / dy ***(x) log \x - y\*%{y) + he . (80) 

Jv Jv 

We first do the integrals in s and s'. Then, the integral in (3 gives a term antisymmetric in x and y, which 
is discarded, and a single symmetric term proportional to the delta function 8(z(x) + z(y)). Using the 
identity 

we find after some algebra 

A3 = — 2m 2 log(m) dx dy log |x — y\ 5(z(x) + z{y)) 
Jv Jv 

2 , t \ f°° r n 1 dxi dxy . . 

= 2m log(m) / dz } log xj — xy — — . (82) 

In the last equation we have written xi(z) as the only solution of z(x) = z in the I th interval, and xi(z) 
as the only solution of z(x) = —z in the I th interval. For a single interval of length L this gives 

A3 = — \og(L) log(m) m 2 L 2 H — log(m) m 2 L 2 , (83) 
3 9 

but for more than one interval A3 does not depend simply on the total length Lt = ^2(bi — ai)- In the case 
of a single interval the entropy can be expressed in terms of a solution of a Painleve ordinary differential 
equation [T3]. The expansion for small mass obtained here coincides in this case with the one (more 
easily) obtained with the help of this differential equation. 

A remarkable feature of the massless fermion entropy (|63p is that the mutual information I(A, B) = 
S{A) + S(B) - S(A U B) is extensive in this case, I (A, B U C) = I (A, B) + I (A, C) [13] . The leading 
correction (|75p for small mass ~ m? log(m) 2 also has extensive mutual information. The first non extensive 
contribution is A3 at order m 2 log(m) |14j . 

4.3 Massive case: long separation distance 

The case of two sets A and B which are separated by a distance r large with respect to their sizes can be 
treated as in [14J. Following [2] the resolvent can be expanded perturbatively, 



R v (f3) = R^(f3) - R^iftdR^iP) + R^ifrdR^iftdR^iP) - ... . (84) 
Here Ry{(3) is the unperturbed resolvent, 

iwi - ( R T z,) ) ■ (85 > 

corresponding to infinite separation distance, and C\ is the field correlator between points in A and B 

TV <«> 
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When r is large C{r) is small and the expansion (I84p makes sense. The correlator between points in 
one set and the other can be considered as a constant in the lowest order approximation (up to correction 
of the order of the cocient of the sises La,b of A and B with |f|). Thus, we have 



The kernel 1 j 4 i b(x, y) = 1 for any x £ A and y £ B and C(r) is given by (|64|) with |x — y| = r. 

If A and -B are small relative to the inverse mass, we can make use of the massless resolvent in (|85p . 
as the zero order term of an expansion in powers of (mLA,B) for Ra(P) and Rb((3). In this case we have 
for the mutual information 

I(A,B) ~ Qd(3 ((3 - 1/2) [tr (c\r)R^((3)C(r)R^((3) + (A «- £?)) - (/? - -/?)] , (88) 

where the bar over the resolvent and the square of the resolvent means sum over the spatial variables, 
Ox = fx dx J x dyO(x, y). Expanding the resolvents in terms of the eigenvectors as above, we arrive at 
the leading term 

I(A,B) ~ im 2 L? Lf (K%(mr) + Kf(mr)) (1 + 0(L A>B /r) + 0{m A L% B ) + ...) . (89) 

In this formula we have to use La,b "C r, m~ 1 , but mr can have any value. Note that (j89H also shows an 
extensive mutual information at the first order in the long distance expansion. 



5 Concluding remarks 

We have shown that the internal dynamics for a multicomponent region and a massless fermion in two 
dimensions is non-local. However, it displays what can be called a quasilocal behavior, since it only mixes 
the fields on a finite number of one dimensional trajectories. The n related trajectories are determined 
by the parameter z which moves at constant velocity. It is interesting to note that the relation —e z = 
Q(x — a)/n( x — can be thought as a "uniformization" transformation [15]. This is the many to 
one transformation which is analytical except at the points a, or bi, and maps each interval into the 
real line. Since this transformation is determined by analyticity properties alone, one could wonder if a 
similar phenomena of quasilocality might also occur for other conformal field theories in two dimensions. 
The factors (|56p correspond to the standard general conformal transformation of the (primary) fermion 
operators. The picture is that of a single variable z with uniform motion, and the n trajectories generated 
by the multivalued transformation. The particularities of each theory might then be reflected in the 
dynamics of the field rotations between the different trajectories. 
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